Non-equilibrium local pair potential enhancement 
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We propose a new mechanism for the local pair potential enhancement with the help of elec- 
tromagnetic control fields. The mechanism is based on the creation of non-equilibrium, spatially 
localized Bogoliubov quasiparticle excitations, which result in a significant enhancement of the local 
pair potential and the local transition temperature T c . 



A properly cooled Fermi atomic gas with an attractive 
interaction between atoms can undergo a superfluid tran- 
sition analogous to the BCS transition in superconduc- 
tors This theoretical prediction motivated many ex- 
periments to study the superfluid properties of cold Fermi 
systems [HHf along with additional theoretical work (see, 
for example, Q and the references herein). Experimen- 
tal studies have been carried out on the manipulation of 
cold Fermi gases [H based on the control of the critical 
transition temperature T c . Control has been achieved 
by application of an external magnetic field 0, [B[ , which 
can modify the effective interaction between atoms due 
to the Feshbach resonance [|| . The enhanced interatomic 
interaction may produce a stronger pairing potential and 
a higher T c . 

In this work we propose a new nonequilibrium mecha- 
nism for pair potential enhancement, which is not based 
on the direct control of the effective interaction between 
fermions, but rather relays on control of the spatial lo- 
calization of the quasiparticle density in the system. The 
local pair potential A(xo,io) at a given position Xq and 
time to is proportional to the product of the local ampli- 
tudes of the quasi-electron and quasi-hole (Bogoliubov) 
excitations v*(xq, to)ui(xo, to) 0- If an external control 
field is chosen in a such way that it drives a significant 
amount of the quasiparticles to be localized at a certain 
moment in the vicinity of xo , this can lead to an enhance- 
ment of the local pair potential A(xq, to)- Therefore, the 
pair potential can be enhanced locally in a target volume 
rather than in the whole system, and this enhancement is 
achieved along with decreasing of the quasparticles den- 
sity and the local pair potential elsewhere. This will in 
effect, result in increase of the local effective critical tran- 
sition temperature T c (xo,io)- A simple estimate of the 
enhancement of the transition temperature at a given 
moment to can be obtained using the BCS expression for 
bulk superconductors T c (xo,to) ~ A(xo, to) /(1.76k b)- 

The suggested mechanism should be contrasted with 
the known Eliashberg mechanism of the pair poten- 
tial enhancement in superconductors Q. Since the 
time-dependent electromagnetic field may create a non- 
equilibrium distribution of quasiparticles in the super- 



conductor, it may lead to unoccupied states at the Fermi 
surface at the gap edge that effectively increases the gap 
value. However, the Eliashberg mechanism does not re- 
sult in a spatially localized enhancement of the quasipar- 
ticle density. The Eliashberg mechanism is rather weak 
effect, resulting in a relative increase of the transition 
temperature T c by an order of 1% , and it was observed 
and extensively studied in the cases of electromagnetic, 
acoustic or tunneling processes 

To illustrate our approach, let us consider the inter- 
action a cold atomic Fermi gas with a time dependent 
control field. A nonequilibrium state of the system is de- 
scribed utilizing time-dependent Bogoliubov-de Gennes 
(TDBdG) equations for an inhomogeneous system Q: 

ih dUn ^ l) =Hu n (r, t) + A(r, t)v n (r,t), (1) 

i% dv n (v,t) = _ HVnM + A *( M K(r,t), 
at 

with H = H + U(r, t) + V(r, t) - fi, H defined by 

Ho = -^ 2 + W(r). (2) 

Here m is electron mass, fi is the chemical potential, 
W(r) is an external trapping potential, and V(r, t) is 
the control field. U (r, t) is the Hartree-like mean field 
local potential, given by 

U(r,t) = -A)^[KM)| 2 /« + Mr,t)| 2 (l - /„)](3) 

n 

And A(r, t) is the quantity of our interest: the local pair 
potential, which is 

A(r,t)=2?oX;u»(r,fK(r,f)(l-2/ B ). (4) 

n 

Here Do is the effective attraction interaction coefficient. 
For example, for isotropic attractive interactions between 
neutral atoms in a trap Do = Airh 2 /ma s , where a s is the 
atomic scattering length. 

To illustrate our approach let us consider a system con- 
sisting of a cooled atomic Fermi gas trapped in a potential 
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VF(r). The trapping potential can be realized, for exam- 
ple, using the induced dipole potential laser trap [l~3j |. 
In this case the potential is proportional to the time av- 
erage of the local laser intensity I(r). The intensity of 
laser can be modulated, which is equivalent to adding 
a time-dependent external control potential V(r, t). For 
simplicity we assumed that the system is elongated and 
that the control field — VV(r, t) is linearly polarized along 
the same direction. Thus, the one dimensional descrip- 
tion is adequate. For simplicity we set the trapping po- 
tential W{x) to be a square well of the length L with 
infinite walls. However, the suggested control strategy 
will work for other types of anharmonic trapping poten- 
tials with non equidistant transition frequencies. The 
system is assumed to be initially in the ground state. 
We used a mesh with m = 16 and m = 32 equidistant 
discretization points for the spatial variable to perform 
numerical integration of Eq.([T]). We test our numerical 
routine by calculating the interaction of a resonant field 
with a single particle in an infinite square well poten- 
tial and compared with the analytical solution for the 
Rabi oscillations between two different states. We found 
an excellent agreement between numerical and analytical 
solutions. 

Optimal control of the local pair potential we formu- 
late as a search for a field V(x,t), which maximizes the 
average value of the pair potential < A > in a target 
area [xq — €q,xq + eo] over a given time interval t$: 

< A >= (2t eo) _1 / / \A(x,t)\dxdt. (5) 

JT—to Jx -e 

In order to understand how one can control the quan- 
tity Eq.©, we do the following analysis. The TDBdG 
equations Eq.([T]) can be approximately solved in the ab- 
sence of the control field (V(x, t) = 0), assuming constant 
pair A(x) = Ao and Hartrcc U(x) — Uo potentials. The 
solution in this case is 

u n (x,t) sa u n ifj n (x)e~ t K , v n (x, t) m v n ip n (x)e~ % K (6) 

with ul = \ (l + = h( 1 ~ where (1 = 

fi — Uq. The eigenenergies are E n = \J (e„ — /2) 2 + A 2 ,, 
where e n and ipn(x) are the eigenenergies and eigen- 
functions for the stationary Schroedinger equation with 
the Hamiltonian Eq.([2]). A control field that maximizes 
< A > should drive at least some of the amplitudes u n , v n 
from their initial state u n (x, 0), v„(x, 0) to a nonequilib- 
rium state at the target time u n (x,T),v n (x,T) (assum- 
ing t <C T), which will have a "bump" at, or in a close 
vicinity of Xq. 

As initial guess for the control field V(x,t) we choose 
one that drives the initial quasiparticle wavefunctions 
with the lowest energy E\, 0), V\{x, 0), toanonequi- 
librium localized state u±(x, T), i>i(x, T) cx g(x), where 
gix) = (x — xq)c~ < -- x ~ Xo,> l a , which has two maximums 



near xq = L/2. The analytical solution for this optimal 
control field in the dipole approximation can be written 
using the results for localization of a particle in an infinite 
well potential obtained in [Io| : V{x,t) = xf(t), 

N 

/(*)= Yl V n cos(u n t/h), (7) 

n=2,4,.. 

where huj n = E n - Ei, V n = a„ = J g(x)ip n (x)dx 

are the expansion coefficients for the target nonequi- 
librium state g(x) in the basis tj) n {x), and d n 
j ipi(x)xip n (x)dx are the corresponding transition dipole 
matrix elements. For symmetric trapping potentials 
W(x) the dipole matrix elements d n = if n is odd, 
therefore we have chosen the target state g(x) to be an 
antisymmetric function with respect to x$. For an in- 
finite well potential e n = 2 mL 2 ' an< ^ ^ e ma t r i x ele- 
ments can be also calculated analytically for even val- 
ues of n : d n = — 4L(cos(7rn)n + n)/(n(l — 2n 2 + n 4 )). 
The expansion coefficients in the limit of small width of 
the target state a <C L, asymptotically approach a n w 
exp(-a 2 (l + n) 2 /4)(l + exp(a 2 n)(n - 1) + n)/4. 
In the present calculations we have chosen Dq = 
10~ 3 Eq. The initial states u n (x,0),v n (x,0) were not 
calculated self-consistently, instead we used the approx- 
imate analytical solution Eq.®. The first six terms 
in Eq.© are used as a starting approximation for the 
control field. We assumed a = 0.1L and set the con- 
trol time interval T = SOnTi, where T x = ^ = 

i is the time period for a particle occupying the 
ground state of an infinite well potential. For this rel- 
atively long control interval the amplitude a quasipar- 
ticle, which occupies the ground state of the potential 
will make about 150 oscillations. Therefore, the Ro- 
tating Wave Approximation, used in the derivation of 
Eq.© in [13 is justified. The optimal control will re- 
sult in relatively slow Rabi oscillations of the occupa- 
tion numbers |a„(£)| 2 ~ |a„| 2 sin 2 (J^), n > 1 [l(|, where 
a n{t) = J ui(x,t)ip n (x)dx. Similar dynamics can be ob- 
served for the quasi-holes b n (t) — J V\(x, t)ip n (x)dx. 

The amplitudes V n and the frequencies w„ of the con- 
trol field were used as an initial input for a black box 
optimization using simulated annealing algorithm. As 
the fitness function for the optimization we used the 
averaged pair potential Eq.([5]) on the spatial interval 
[L/2 - e ,L/2 + e Q ], with e = 0.1L, and to = CUT. 
The simulated annealing optimization resulted in addi- 
tional « 15% improvement of the averaged magnitude 
of the pair potential < A > with respect to the en- 
hancement, obtained using the non-optimized analytical 
solution Eq.([7]). To reduce the amount of the meshing 
points we set the renormalized chemical potential jl to 
zero, therefore limiting the amount of the nodes in the 
initial quasiparticle amplitudes. We considered the evo- 
lution of maximum 6 (3 quasi-electron, and 3 quasi-hole) 
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quasi particles. Such a small number was chosen be- 
cause the characteristic single particle frequencies in a 
square well potential grow quickly (as n 2 ) that makes nu- 
merical integration of TDBdG equations extremely slow. 
The resulting nonequilibrium dynamics of the pair po- 
tential |A(x, i)| is shown in FigjTJ Note we assumed 
zero temperature. In Fig. 1(a) one can see the gradual 
spatial localization of the pair potential near the cen- 
ter of the infinite well potential xq = L/2, which occurs 
at the end of the control interval. The local enhance- 
ment of the pair potential and the transition temperature 
T c {xq, T) oc A(xq, T) is about 300% compared to the ini- 
tial state. In FigQJa) one can see that the enhancement 
is achieved at the price of decrease of the |A(x, t)| out- 
side the target area, closer to the walls of the trapping 
potential. In Fig.l (b,c) we have shown similar simula- 
tions for optimal control of four and six quasiparticles, 
and the resulting dynamics of the pair potential. As 
in the case of just two quasiparticles, the pair potential 
shows a strong enhancement near the center of the well. 
Note, the control field V(x,t) also drives other quasi- 
particle amplitudes u n (x,T), v n (x,T) with n > 1, but 
this does not change much the localization picture. One 
has to note that the enhancement is relatively weaker in 
the latter two cases shown in Fig.l (b,c). The approx- 
imate analytical solution for the amplitudes a n (t) and 
b n (t) has the period of AT, therefore, for the controlled 
system |a„(i)| 2 and |fr n (t)| 2 will return to their initial 
states at t = 2T. In real systems this return will never 
be perfect, and it will correspond to the loss of coherence 
in the driven system and the energy dissipation. For 
strongly interacting quasiparticles the driven dynamics 
can significantly depart from the approximate analytical 
solution. In superconducting materials the pair poten- 
tial A is relatively small compared to the Fermi energy, 
and can be treated as a perturbation. However, this is 
not always the case, for example, in cold atomic Fermi 
gases. Relatively strong coupling between functions u n 
and v n in the presence of oscillatory field may lead to 
chaotic dynamics. In this case numerically stable control 
techniques of chaos are necessary [ill ] . 

The proposed mechanism may be applied not only to 
cold atomic Fermi gases, but also in solid state quasi- 
one dimensional superconducting systems. However, the 
thickness of the system should not be too small, since 
the gap enhancement can be suppressed by the phase slip 
phenomena f]~2j ] . which is not included in our mean field 
description. One may consider an experiment with quasi 
one dimensional superconductors, such as a nanowire or 
a bundle of doped carbon nanotubes 14]. One may de- 



sign an optimal control experiment, when the quantity 
of interest will be non-equilibrium conductivity of the 
bundle at a temperature close to T c . Making analogy 
with the Eliashberg effect [ill ], one may expect that a 
nonequilibrium perturbation of the system may increase 
the effective critical temperature (as well as the instanta- 
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FIG. 1. Spatial localization and enhancement of the pair po- 
tential |A(a;,£)| (arbitrary units) using an optimized exter- 
nal control field. The optimization was performed with the 
target of maximizing the fitness function in Eq.©. (a) two 
quasi-particles in the infinite well potential (described by the 
amplitudes u\(x,t) and Vi(x, t)) (b) four quasi-particles, (c) 
six quasiparticles. 



neous pairing potential). Using an optimal control field 
one may drive periodically the quasiparticle density to be 
maximal along the bundle near its center, that will re- 
sult in a considerable drop of the resistance of the bundle 
at times of the spatial localization. In order to achieve 
this, the control field should have its polarization, which 
is transverse to the bundle. 

The control field Eq.([7]) can be used in the experi- 
ment as an initial guess. Then one may use the mea- 
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sured nonequilibrium conductivity of the bundle at the 
target moment T as a fitness function, in a similar fash- 
ion we used Eq. ([5]) for the simulated annealing optimiza- 
tion. Let us consider a bundle of boron or alkali-atoms 
doped carbon nanotubes of width L = 500nm [14j. We 
assume the level of doping such that the carriers are hav- 
ing the wavelength of the order of 4 nm. As the ini- 
tial guess for the control field we may choose a field that 
drives the initial quasiparticle wavefunction with the low- 
est energy u n (x) — W250 oc sin(2507rx/L) to the target 
U25o(x,T) oc u n (x,0)g(x), where the function g(x) is de- 
fined above. Note the control will be performed rather 
over the envelope of the quasiparticle wavefunctions. At 
the time T the target wavefunction can be represented 
as u 2 5o(x,T) = Y,j=i,a,.. X) n =250±j a n u n (x,0). The ex- 
pansion coefficients of the target wavepacket a n can be 
used in the solution Eq.([7]). The same control is applica- 
ble for the quasi- hole wavefunction V25o(x). Even for a 
relatively modest density localization with a = 0.3L for 
just one pair of quasiparticles with the lowest energy, it 
will correspond to s» 300% increase of the term w^so^aso 
and the overall local increase of the pairing potential by 
« 1%. One can use the adaptive optimization approach 
to further improve the localization of quasiparticles [l6| . 
There may be some possible obstacles for the proposed 
experiment. First, relatively strong control field may re- 
sult in strong energy dissipation. At the same time the 
control field can not be too weak, since it will increase the 
duration of the optimal control interval T. The duration 
of the control interval should not be too long, because the 
loss of coherence, for example, due to collisions of cold 
atoms, may reduce the efficiency of the proposed scheme 
17|. 

To conclude, we presented a new approach to con- 
trol the local enhancement of the pair potential. The 
enhancement is achieved through a nonequilibrium spa- 
tial localization of the quasiparticle density in the sys- 
tem. This mechanism contrasts with the spatially ho- 
mogeneous Eliashberg mechanism, which usually em- 
ploys a field with only a single frequency. Our method 
is based on coherent control, which requires a multi- 
frequency field with independently tuned field ampli- 
tudes. The suggested control scheme will be effective 
on the time scale, shorter than the effective decoher- 
ence times, which for Fermi gases can be of the order 
of milliseconds. The time-dependent Eq.([T]) is valid for 
weakly nonequilibrium conditions, since A(r, t) is defined 
by using the equilibrium values of the Fermi distribution 
/„. Most significant, the electromagnetic field should not 
have components with the frequencies over the double 
gap size, huj max > 2|A|, because high energy photons 
can easily break Cooper pairs. Assuming the equilib- 



rium pairing potential A « 0.0015eV in the above de- 
scribed example with the bundle of doped carbon nan- 
otubes, one may estimate that the transition frequen- 
cies between the nearest levels at the Fermi surface (as- 
suming the infinite well potential model) will be about 
AE « 0.0007eV. Therefore operating with a radio fre- 
quency field with fko max < 3 x 10 _3 eE one still will be 
able to create a linear superposition of states U25o(x, T) — 
a25i"25i (x, 0) + 0253^253(2;, 0) -I- 0255^255(0;, 0), which has 
a considerable "bump" near the center of the bundle. We 
would like to emphasize that direct microscopic simula- 
tions using Eq.lU) will need at least a thousand spatial 
discretization points that makes the problem currently 
computationally intractable. 
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